Gravitational-wave bursts from soft gamma-ray repeaters: Can they be detected? 
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In this letter we suggest a scenario for simultaneous emission of gravitational-wave and 7-ray bursts 
(GRBs) from soft gamma-ray repeaters (SGRs). we argue that both of the radiations can be gener- 
ated by a super-Eddington accreting neutron stars in X-ray binaries. In this model a supercritical 
accretion transient takes back onto the remnant star the disk leftover by the hydrodynamic instabil- 
ity phase of a low magnetized, rapidly rotating neutron star in a X-ray binary system. We estimate 
the rise timescale At c — 0.21 ms, minimum mass accretion rate needed to trigger the 7-ray emis- 
sion, M\ — 4.5 x 10 28 g, and its effective associated temperature T e ff = 740 keV , and the timescale 
for repeating a burst of 7-rays Atr = 11.3 yr. Altogether, we find the associated GW amplitude 
and frequency to be h c — 2.7 x 10~ 23 /(Hz) 1 ^ 2 and f gw — 966 Hz, for a source distance ~ 55 kpc. 
Detectability of the pulses by the forthcoming GW anntenas is discussed and found likely. 

PACS numbers: 04.30Db, 97.60. Jd, 98.70Rz 



With the advent of the gravitational-wave astronomy 
a new window in physics is about to be opened. The 
first GW signal, whenever detected with the new gener- 
ation of observatories such as LIGO, VIRGO, GEO600, 
TAMA 0| , and / or the fourth generation of resonant- mass 



detectors like TIGAs [2 
modern astrophysics [7 



will mark a breakthrough in 
Most theoretical studies to 
address the prospects have been concentrated in sources 
like coalescing NS binaries HfL0| and rapidly rotating 
NSs 1 11 I^|l4|,|l5|] . Statistical estimates of event rates, 
for both sources, point to an occurrence rate ~ 3/yr, 
in the best cases, for distances up to the Virgo cluster 
pG| p"8[ . On the other hand, neutron star hydrodynam- 
ical instabilities, like that one studied by Houser, Cen- 
trella & Smith (1994) (HCS'94) and Lai & Shapiro (1994) 
jl2],[l4| , should occur quite frequently in nature and we fo- 
cus on the post-hydrodynamical instability phase in such 
systems, with the perspective of modeling SGRs, which 
may also be interesting GW sources Ufl . 

Because of the physics underlying the exact nature of 
soft gamma-ray repeaters (SGRs) is still an unsolved puz- 
zle in high energy astrophysics, the present study is a 
novel ingredient in our search for understanding these ob- 
jects. Together with, it is the realization that the foresee- 
able detection of gravitational waves (GWs), as discussed 
above, may help underpinning various open problems in 
both lines of research, yet. The possible confluence of 
these widely different trends in modern astronomy is an 
interesting perspective. Motivated by such a worthy pos- 
siblity, we propose that a neutron star accretes matter 
from a companion in a X-ray binary system until becom- 
ing a millisecond spinning star. Then, centrifugal effects 
drive equatorial mass ejection, and form a thick disk sur- 
rounding the NS. Various instabilities in this dense struc- 



ture forces it to fragment into clumps that are later on 
capture back onto the star. Our model allows us to es- 
timate basic important features of SGRs as its risetime, 
minimum mass accretion rate needed to trigger the 7-ray 
emission and its effective associated temperature. The 
time interval to becoming a repeating 7-ray burster is a 
worthy by-product of our scenario, as shown below our 
results are in quite well agreement with observations. 

Since most of SGRs observed fluxes indicate a luminos- 
ity ~ iq4i-42 er gj s ( or probably higher, ~ 10 44 erg/s 
p9[); well above the Eddington luminosity, it is nat- 
ural to consider these binary systems in a supercriti- 
cal accretion regime. In our SGRs model we envision 
a disk-like structure as being a donut-type thick encir- 
cling disc with characteristic dimensions, (see Table-1): 
Aw = 2 x R ns as the disk width; the difference between 
the external radius Rd', assumed to be that one derived 
for the final structure in HCS'94 simulations, Rd ~ 5R ns , 
and the internal radius, supposed to be the tidal radius, 
Rude ~ SRns- We also assume that the disk density is 
near that one for the HCS'94 remnant disk-like structure, 
i. e., pa ~ 5.7 x 10 12 g/cm 3 p2| . For parameters of our 
NS model see Table-2 |l§4 

Fluid structures, like this accretion disk-type object, 
may undergo various instabilities, depending on the dy- 
namical conditions (Jeans, Rayleigh- Taylor instability, 
viscous). Provided the system is driven by Jean's in- 
stability, clumps of matter around the central compact 
body are produced. The lengthscale of these clumps fol- 
lows from the Jeans' wavenumber 
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where Kj is the wavenumber of the perturbation in 



1 



TABLE I. 

Parameters of our disc-like structure: mass Md, density 
Pd, width Aw, external radius Rd, tidal radius Rude, or- 



biting radius R or b, 
dius R acc . 



disc height-scale AH and Afven ra- 



Parameter 



Value 



M d [M Q ] 
Pd [g/cm 3 ] 
Aui [cm] 
Rd [cm] 
Rtide [cm] 
Rorb [cm] 
AH [cm] 

Race [Rq] 



4.0 x 1CT^ 
5.7 x 10 12 
2.82 x 10 6 
7.05 x 10 6 
4.21 x 10 6 
5.64 x 10 6 
5.0 x 10 4 
1.0 x 10 3 



the matter, ad — (-^j) 2 is the sound speed in the struc- 
ture and pd its mass density. Hence, assuming that ad ~ 
10 _1 (c) cm/s, we find Kj ~ 7.3 x 10 -7 cm -1 , and there- 
fore determine an associated length, Xj ~ 8.7 x 10 6 cm. 
We attribute the production of a burst to the fall of a set 
of clumps back onto the NS (see below) . 

In this scenario the rise timescale and the time-delay 
between peaks in the soft 7-ray bursts spectra are deter- 
mined as follows: when fragments fall back onto the NS 
surface, with balistic trajectories, each blob is stretched 
to a needle shape 20 1 . The time interval for this collision 
(blob-NS) to occur would be 
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here M ns and R ns are the mass and radius of the 
NS, respectively. Therefore, we argue this time interval, 



At r 



2.1 x 10 s, is that one characteristic of rise- 



times in SGRs. For comparison, it is worth noticing that 
the 790305 b event had a rise time < 0.25 ms @,|(|. 
The characteristic time between intervening peaks in a 
given burst structure, (the timebreak), could roughly 



be inferred from r s 



= K-j x V rb and turns out to 



be 



168.2 ps. Here V t 



orb 



(0.273 c) 



cm/s, is the orbital velocity of each clump before being 
accreted. This value is over the GRO temporal resolu- 
tion: ~ 100 ps H^Jl^], and therefore may be observed. 
In clumped structures such as these strong turbulence 
effects should play an important role in driving the disk 
hydrodynamics. We may take these effects into consid- 
eration in our present simplified description throughout 
the parameter /3 QCC , which tells us about the efficiency of 
viscosity in governing the structure. 

After losing enough energy and part of its angular mo- 
mentum through those mechanisms, the clumped disc- 
like mass distribution falls onto the NS on a free-fall 
timescale. When the blob hits on the NS surface it trig- 
gers a strong thermal-like burst of 7-rays. The charac- 
teristic temperature of the burst can be estimated from 



equating the luminosity of accretion to the radiation 
power. This yields the equation 
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M\ « ^npd{XAwAH), with pd being the disc den- 
sity; T pea k is the radiation temperature, and M remn and 
Rremn represent the mass and radius of the remnant com- 
pact star. We define M\ ~ being Tff = (Gp)^ 1 ^ 2 ~ 
548ms. The equivalent accretion rate per year yields 
M\ ~ 1.1 x 10 5 M^/yr. Assuming that the hard com- 
ponent of GRB 790305 b is that one characteristic of 
all SGRs, it can be seen that most of SGRs observed 
thusfar would have had associated mass accretion rates, 
M\ ~ 4.5 x 10 28 g/s. We shall assume a value of the pa- 
rameter (3 aC c ~ 0.70 for the case of accretion onto NSs.f 

TABLE II. 

Our NS model parameters: period P, mass M ns , radius 
R ns , angular momentum J, rotational to kinetic energy 
ratio f3. For initial we mean at the advent of the insta- 
bility, and for final the stable stage after it. 



initial 



final 



P [ms] 
M n3 [M Q ] 
R„s [cm] 
J [gem 2 /s] 
P 



0.900 
1.47 
1.37 x 10 6 
9.84 x 10 

0.2738 



,4!) 



0.98 

1.42 
1.41 x 10 6 
8.78 x 10 

0.27 



,411 



By using parameters derived from our scenario (see 
Table 1), we obtain for the burst peak temperature a 
value Tpeafc ~ 2.6 x 10 10 K, or equivalently, T pea k ~ 
2.2 MeV. 

Observations of SGRs J22| have shown that in GRB 
790305 b there was a soft spectral component which con- 
tained most of the energy and 90% of the photons be- 
tween 30 and 2000 keV. Nevertheless, for blackbody 
models of the emission processes the effective temper- 
ature, i. e., that one is measured by the detectors, 
Teff ~ Tpeak/3. Thus, we infer from our model the value 
740 keV which is comparable to the hard component ob- 
served in SGR 790305 b, 520 ± 100 keV [|fj). This leads 
us to conclude that the scales of temperatures predicted 
and observed agree with each other. Higher accretion 
rate would yield in stronger 7-ray bursts. That it would 
be the case, for example, for the expected accretion rate, 
- 8.9 x 10 30 g/s, in the HCS'94 ejected disk-like matter; 
we expect the maximum of the emission in that situation 
to peak at ~ 3.8 MeV and the associated GW pulse in 
that event to be slightly stronger than the one computed 
for SGRs in the present estimate (see below). 

Viscosity torques remove angular momentum forcing 
the disk eventually to fall down onto the star after a 
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formation. The compact body recovers angular momen- 
tum, and rotates more quickly. The amount of angular 
momentum restored, after the accretion of all the clumps 
of matter, reads 



Pace x N X M A X V orb X R orb ~ A J r 



(0.4) 



where N represents the number of clumps, ~ Kj x 
2nR or b ~ 5, R rb defines the radius for bound orbiting, 
~ 4.0i? Ils , (HCS'94), and AJ rec the recovered angular 
momentum. With our assumed parameters for the ac- 
creted material we get AJ rec ~ 5.2 x 10 48 gcmPs -1 . 

The timescale to repeat a 7-ray burst can be obtained 
from the amount of the angular momentum to be ac- 
creted from the main accretion disk (the donor star) in 
order to reentering into a new unstable phase. To esti- 
mate such contribution we first calculate the value of the 
energy parameter, f3, after the transient. This yields 
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= 0.27, 
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where we have introduced the factor / to measure 
the departure of the star from the dynamical instability, 
which lies the interval 0.3 < / < 0.6, and parametrizes 
the actual star's non-sphericity and polytropic index of 
its equation of state in the post-hydrodynamic instability 
phase. Here the subscript word stable stands for secular 
phase. The lower limit for / indicates the star has en- 
countered the secular instability. The higher one guaran- 
tees the achieving of the dynamical limit. In the present 
situation we expect this factor / to be ~ 0.6. The /3 value 
for the stable phase, in eq.(5), determines the amount of 
angular momentum to be gained from the primary star 
at the characteristic rate of accretion in such systems 
M ~ (10~ 7 - 10~ 5 ) M Q /yr |IJ. 

The next step is to determine Atr, the time to re- 
peat a burst. Within a burster phase, we use the rela- 
tion between the inward angular momentum, J + , and 
the amount of angular momentum to be injected from 
the main accretion disk, A Jr (so that the star becomes 
unstable again at ~ 0.2738), to write 



J + x At r ~ (3 acc x AJ R , 



(0.6) 



where Atr defines the timescale for bursting again, 



and J" 1 



M d (GM stabl 

e-^acc) 
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21|. Here Md defines 



the rate of mass accretion from the main disk, R acc cor- 
responds to the distance the accreted matter will come 
from, the Afven radius for the low B NS in the system, 
~ lO 3 i?0. Thus, if we impose for our model NS an ac- 
cretion rate Md ~ 10~ 6 M@/yr we find Atr ~ 11.3 yr. 
This result is in the observed ballpark figure of SGRs ac- 
tive phase timescales ~ 10 years. Particularly, the source 
SGR 1900+14 was observed for the first time in 1979 
|§. The SGR 1806-20 was observed in 1987 Q. Both 
became active again in the nineties as confirmed by ob- 
servations from BATSE (CGRO) and ASCA @Mp3. 



SGR 0526-66 was detected in December 1981 and April 
1983, according to observations from the spacecrafts Ven- 
era 13 and 14 [^(J. In this picture, the accretion rate 
must become strongly supercritical in the last source, 
M ~ 10~ 4 M Q /yr, so that the time scale would turn 
out to be ~ 0.11 yr, close to observed time intervals of 
gamma-ray recurrences in some SGRs, which are ~ 39 
days [Q. This result is along the lines of King et al. 
(1997) model [ p5| , and suggests that it can find SGRs 
systems undergoing episodes of irradiation-driven super- 
critical mass transfer (IDSMT) during short time spans, 
depending on binary orbital parameters, with recurrences 
occurring at time steps as short as ~ 40 days. 

Now, for the case of the theoretical HCS'94 NS rem- 
nant, and using the same distance and rate of accretion 
as above, with a stable mass for the NS ~ 1.38M Q , we 
obtain for Atr ~ 2.4 x 10 3 yr. Hence, if the accre- 
tion process were supercritical, as in systems undergo- 
ing IDSMT, the resulting timescale for repeating a burst 
would be ~ 24 yr, what is nearly two orders of magnitude 
longer than that one for Soft Gamma-Ray Repeater Low 
Mass X-Ray Binaries (SGR-LMXBs), just estimated in 
the last paragraph. 

Finally, we wish to determine the characteristic GW 
amplitude associated with the burst of 7-rays. In our 
model, it is the sudden accretion of matter what trig- 
gers the emission of a pulse of GWs during the transition 
phase to a stable rotating NS, through the fluid modes /. 
To estimate the GW burst amplitude we can compare the 
GW flux received at Earth to the total gravitational lu- 
minosity radiated away by the source (§] , extracted from 
the gravitational potential energy of the orbiting disk-like 
structure 
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where AEqw represents the energy flowing into GWs, 
and is defined as AEqw — ALgw x At c , with At c de- 
fined in eq.(2). The gravitational luminosity ALqw is 
given by M 



AL GW = K x 



GM\ L 



2 r 4 
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where K is a constant: 10~ 2 < K < 10 |§7). We 
assume the value K = 0.5, the mass being M\ = 
3.6 x 10 30 g; the one that triggers the emission of the 
highest energy photons, and the parameter L is the radius 
of the orbit of that mass around the remnant, ~ 4i?„ s . 
Introducing these figures in the equation 



h c = 



AG / 9 
^(AWx(At c ) 2 



-1 1/2 



we get for the GW amplitude the value h c 
10~ 22 , for a source distance, D, of 55 kpc. 
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In order to be detected by the TIGAs network, the 
SGRs characteristic GW amplitude, when normalized 
to its frequency Fourier transform (the square root 
of the characteristic spectral frequency distribution), 
our estimated GW-amplitude reduces to h c — 2.7 x 
10" 23 /(iJz) 1 ^ 2 . We propose that the GW-frequency 
should be correlated with the timescale during which half 
of the infalling clumps of matter impact the NS surface, i. 
e -, fgw ~ 1/(5 x At c ) = 966 Hz. This timescale specifies 
the elapsed time until a " clean" GW signal from the sys- 
tem is settled down in the GW detector. The frequency 
sensitivity bandwidth being Af gw = 2fj,f gw ~ 30 Hz, 
with [i the inverse of the two-mode transducer mass ratio. 
Such frequency bandwidth is expected to be achieved for 
most of the planned resonant-mass detectors |||). The 
frequency-averaged GW amplitude is certainly under the 
quoted TIGA's sensitivity cut off (see Figure 1). 

As a check of the calculation just done, we can com- 
pute the GW characteristic amplitude with the Wag- 
oner's (1984) equation 11 rescaled for the actual accre- 
tion mass rate, M\, in our scenario 
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This yields h c = 3.0 x I0~ 23 / \H z) 1/2 '. The calcula- 
tion for the whole range of frequencies (1-3000) Hz are 
presented in Figure- 1. All of the plotted values were nor- 
malized for the frequency spectral distribution (Hz) -1 / 2 
of the outgoing GW burst, in order to be compared to 
the TIGAs network and LIGO sensitivities, as shown in 
Figure- 1. 

Our model is based on the hypothesis that there ex- 
ists a peculiar class of LMXBs where higher accretion 
rates are at work. The burst duration in SGR-LMXBs 
should be associated with the total accreted mass for the 
transient. As shown above, the time span for re-bursting 
can be inferred to from the burst temperature. Dur- 
ing quiescent stages, assuming an Eddington's accretion 
rate and a source distance ~ 55 kpc, their X-ray fluxes 
must be around 3.4 x 10~ 16 erg/(cm 2 s). These fluxes 
are marginally detectable by Beppo-SAX or ASCA X- 
ray detectors. For galactic distances, ~ 10 kpc, those 
sources would certainly be detectable by such capabil- 
ities. For comparison, observed values for the SGRs 
X-ray luminosities are L x = 7 x 10 35 erg/s and L x — 
3 x 10 34 (£>/8fcpc) erg/s, for SGR 1900+14 and 1806-20, 
respectively pq| . 

As a point aside, there are left two further pathways 
for the thick disk to evolve. Firstly, it may not fragment 
apart anyway. Thus, the whole disk infalls at once onto 
the NS triggering, therefore, more powerful bursts of both 
of the radiations. Secondly, had fragments it produces 
very small lenghtscales, then, instead of being observed 
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FIG. 1. 

Locus of the dimensionless GW characteristic peak am- 
plitude, h c (log), plotted against the wave frequency, 
few (linear), for the three up-to-date known soft 7-ray 
repeaters (distance^pc) H~fl ux [erg/(cm 2 s)\ 

): SGR 0526-66 

(55,- 1.5 x 10~ 3 ), SGR 1900+14 (14 - 4.0 x 10~ 4 ) and 
SGR 1806-20 (17,- 2 x 10" 4 ). Also data are plotted for a 
hypothetical strong 7-ray source, 7-flux — 5 x 10~ 2 (top- 
line), placed at 8 kpc. The symbol X represents the GW 
characteristics computed with eq.(9). 



as a 7-rays train of pulses it should appear as persistent 
bursts of hard or even soft X-rays, with characteristic 
timescale — 1 month, and time recurrences far too short, 
— 3 mins, alike to the ones observed by BATSE from the 
source GRO J1744-28 Jj9 30 1 . Sources such as these may 
prove interesting astrophysical labs to test this scenario. 
Also do numerical simulations of thick dense disk neu- 
tron star interactions taking into account relativistic ef- 
fects, or instead observations of accretion X-ray systems 
in which accretion rates — 10~ 6 Mi^fyr are seemingly 
common!, or even such systems as those simulated by 
Nomoto (1986) and Saio & Nomoto (1998) 
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Our main result is displayed in Figure- 1 which shows 
that the "xilophone" of TIGAs, with a bandwidth — 
30 Hz, may observe GWs through the NS fluid modes 
from SGRs, whenever they radiate at frequencies lower 
than 1.4 kHz, for distances as far as the LMC. As an ex- 
ample of the detectability of these signals by ressonant- 
mass detectors it is plotted the pulse computed with 
eq.(7), marked with an (X), which corresponds to h c = 
2.7x 10- 23 /(Hz) 1/2 and a f gw = 966 Hz. Interferometric 
GW detectors like LIGO, for example, may observe such 
systems provided they emit at mid frequencies, < 640 
Hz (first LIGO), and up to 2000 Hz (advanced LIGO), 
for the same distances. Given the potential implication 
for the physics of both SGRs and GWs, detailed studies 
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of the model would be needed. 
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f Note that this choice guarantees us that the dynam- 
ically unstable NS gets rid of sufficient angular momen- 
tum so as to becoming itself secularly unstable (see also 
Shapiro & Teukolsky 1983). 

% More details on tables 1 and 2, and also on suggested 
modelling and observational work will be given in a work 
in preparation. 
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